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The evolution of the magnetic moment and various features of the electronic structure of fee Gd are 
followed to reduced volume V/Vo = 0.125 using the LDA+U correlated band method. The stability 
of the moment is substantial; crude estimates of this signature of a possible "Mott transition" in the 
4/ system suggest a critical pressure Pc ~ 500 GPa. The 4/ occupation is found to increase under 
pressure due to broadening and lowering of the minority states. This trend is consistent with the 
interpretation of x-ray spectra of Maddox et al. across the volume collapse transition at 59 GPa, 
and tends to support their suggestion that the delocalization of the 4/ states in Gd differs from the 
original abrupt picture, being instead a process that occurs over an extended range of pressure. 



I. INTRODUCTION 

The behavior of the 4/ rare earth metals and their 
compounds under pressure has been discussed for 
decades, with the volume collapse transition under 
pressure attracting a great deal of attention. It has 
been known for some time that there are volume col- 
lapse transitions in Ce (15% at 0.7 GPa), Pr (10% 
at 20 GPa), Gd (5% at 59 GPa) and Dy (6% at 
73 GPa), while no significant volume collapses have 
been detected in Nd, Pm, and Sm. The equation of 
state of these metals, and references to the original 
work, has been collected by McMahan et al} High 
temperature experiments^ have seen signatures that 
are likely related to the localized— s-itinerant transi- 
tion, at 50 GPa in Nd and 70 GPa in Sm. 

The question can be stated more generally as: 
what form does the localized — > itinerant transition 
of the 4/ states take, and what is the correct de- 
scription? This transition is intimately related to 
the question of behavior of magnetic moments (A al- 
though the questions are not the same. There have 
been two main viewpoints on the volume collapse 
transition. One is the "Mott transition of the 4/ sys- 
tem" elaborated by JohanssoufS, in which the crucial 
ingredient is the change from localized (nonbonding) 
to more extended states (participating in bonding), 
with an accompanying drop in magnetic tendency. 
The other is the "Kondo volume collapse" view in- 
troduced by Allen and Martin'* and Lavagna et al.^, 
in which the main feature is the loss of Kondo screen- 
ing of the local moment, with a decrease in localiza- 
tion of the 4/ state not being an essential feature. 

At ambient conditions the 4/ electrons form a 
strongly localized /" configuration that is well char- 
acterized by Hund's rules. Under reduction of vol- 
ume, several things might be anticipated to happen. 
At some point the 4/ system begins to respond to 



the non-spherical environment. Initially, perhaps, it 
is just a matter of crystal field splitting becoming 
larger. Then the 4/ orbitals actually begin to be- 
come involved in the electronic structure, by over- 
lapping orbitals of neighboring atoms. The conse- 
quences of this are possible participation in bond- 
ing, and that the orbital moment becomes less well 
defined (the beginning of quenching i.e. the loss 
of Hund's second rule, which has already occurred 
in magnetic 3d systems) . Additionally, the 4/ levels 
can shift and increase their interaction with the itin- 
erant conduction (c) bands (Kondo-like coupling), 
which can change the many-body behavior of the 
coupled 4/ — c system. At some point the kinetic 
energy increase, characterized by the 4/ bandwidth 
Wf, compared to the on-site interaction energy Uf 
reaches a point where the spin moment begins to 
decrease. Finally, at small enough volume (large 
enough Wf) the 4/ states simply form nonmagnetic 
conduction bands. 

Just how these various changes occur, and in what 
order and at what volume reduction, is being ad- 
dressed in more detail by recent high pressure exper- 
iments. Here we revisit the case of Gd, whose vol- 
ume collapse was reported by Hua et al.^ and equa- 
tion of state by Akella et alJ. The deviation from the 
series of close-packed structures below Pc=59 GPa 
and the lower symmetry bcm (body-centered mono- 
clinic) high pressure structure signaled the expected 
onset of /-electron participation in the bonding, and 
Hua et al. seemed to expect that the moment re- 
duction and delocalization of the 4/ states would 
accompany this collapse. 

New information has been reported by Maddox et 
al.^ who have monitored the resonant inelastic x-ray 
scattering and x-ray emission spectra of Gd through 
Pc and up to 113 GPa. They find that there is no 
detectable reduction in the magnetic moment at the 
volume collapse transition, so the volume-collapse is 



only a part of a more complex and more extended 
delocalization process of the 4/ states. Maddox et 
al emphasize the Kondo-volumc-collaps o?'^'^ aspects 
of the transition at Pc. 

The treatment of the 4/ shell, and particularly the 
volume collapse and other phenomena that may arise 
(see above), comprises a correlated-electron prob- 
lem for theorists. Indeed there has been progress 
in treating this volume-collapse, moment-reduction 
problem in the past few years. The issue of the 
(in)stability of the local moment seems to involve 
primarily the local physics, involving the treatment 
of the hybridization with the conduction bands and 
interatomic f — f interaction, with Kondo screen- 
ing of the moment being the subsequent step. Dy- 
namical mean field studies of the full multiband 
system have been carried out for Ceiii*iS and for 
Pr and Nd.^^ These calculations were based on 
a well-defined free-energy functional and included 
the conduction bands as well as the correlated 4/ 
bands. One simplification was that only an orbital- 
independent Coulomb interaction U was treated, 
leaving the full orbital-dependent interaction (fully 
anisotropic Hund's rules) for the future. Density 
functional based correlated band theories have also 
been applied (at zero temperature). Self- interaction 
corrected local density approximation (LDA) was 
applied to Ce, obtaining a volume collapse compara- 
ble to the observed oneii^ Four correlated band the- 
ories have been applied-'-^ to the antiferromagnetic 
insulator MnO. Although their predictions for criti- 
cal pressures and amount of volume collapse differed, 
all obtained as an S=5/2 to S=l/2 moment collapse 
rather than a collapse to a nonmagnetic phase. 

Clearly there remain fundamental questions about 
how the magnetic moment in a multielectron atom 
disintegrates as the volume is reduced: catastroph- 
ically, to an unpolarized state, or sequentially, 
through individual spin fiips or orbital-selective de- 
localization. If the latter, the total (spin -I- orbital) 
moment could actually increase initially in Gd. If 
the occupation change is toward the decrease in 
spin moment (from S—7/2 to S=3) could be more 
than compensated by an orbital moment (L=3). If 
the change is toward the onset of an L=3 orbital 
configuration could oppose the S=3 spin (Hund's 
third rule), leaving a non-magnetic J=0 ion (as in 
Eu^+) even though the 4/ orbitals are still local- 
ized. Still another scenario would be that the in- 
creasing crystal field quenches the orbital moment 
(as in transition metals) and the remaining problem 
involves only the spins. 

Our objective here is to look more closely at the 
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stability of the Gd atomic moment, in the gen- 
eral context of the localized— >itinerant transition of 
the 4/ system under pressure. Consideration of 
the changes in electronic structure under pressure 
go back at least to the broad study of Johansson 
and Rosengren^^ but most have not considered the 
magnetism in detail. We apply the LDA-I-U (local 
density approximation plus Hubbard U) method to 
study the evolution of the electronic structure and 
magnetism as the volume is reduced. Although this 
correlated band method neglects fluctuations and 
the dynamical interaction with the conduction elec- 
trons, it does treat the full multiorbital system in 
the midst of itinerant conduction bands. The re- 
sulting moment vs. volume surely provides only an 
upper limit to the pressure where the moment de- 
creases rapidly. However, we can invoke studies of 
the insulator-to-metal transition in multiband Hub- 
bard models to provide a more realistic guideline on 
when the localized— >itinerant (or at least the reduc- 
tion in moment transition within the 4/ system may 
be expected to occur. The results suggest stability 
of the moment to roughly the 500 GPa region. 



II. ELECTRONIC STRUCTURE METHODS 

In the paper we apply the full potential local or- 
bital codaM (FPLO5.0G-18) to Gd from ambient 
pressure to very high pressure (a few TPa). We use 
the fee structure with space group Fm3m (#225) 
and ambient pressure atomic volume (corresponding 
to the fee lattice constant ao=5.097A). The basis set 
is (core)::(4ci4/5s5p)/6s6p5(i-|-. We use 48^ k point 
mesh and Perdew and Wang's PW92 functionalii for 
exchange and correlation. We have tried both 5.0 
and 6.0 for the confining potential exponent, with 
very similar results, so only the results using expo- 
nent=5 will be presented here. We perform both 
LDA and LDA-I-U calculations (see below). Due to 
the extreme reduction in volume that we explore, 
any band structure method might encounter difficul- 
ties. For this reason we have compared the FPLO 
results on many occasions with parallel calculations 
with the full potential linearized augmented plane 
wave method WIEN2k:^- The results compared very 
well down to V/Vo=0.5, beyond which the WIEN2k 
code became more difficult to apply. We use the no- 
tation V = V/Vo for the specific volume throughout 
the paper. 

We assume ferromagnetic ordering in all calcula- 
tions. The Curie temperature has been measured 
only to 6 GPaji^iSSiSi where it has dropped from 



293 K (P=0) to around 210 K. Linear extrapola- 
tion suggests the Curie temperature will drop to zero 
somewhat below 20 GPa. However, as the 4/ bands 
broaden at reduced volume the physics will change 
substantially, from RKKY coupling at ambient pres- 
sure finally to band magnetism at very high pressure. 
Antiferromagnetic (AFM) ordering does not affect 
the 4/ bandwidth^^ until / — / overlap becomes ap- 
preciable. AFM ordering might affect some of the 
quantities that we look at in this study at very high 
pressure, but such effects lie beyond the scope of our 
present intentions. 



A. LDA+U Method 

For the strength of the 4/ interaction we have used 
the volume dependent U{V) calculated by McMa- 
han et al.^ which is shown below. Due to the lo- 
calized 4/ orbital and the large atomic moment, we 
use the "fully localized limit" version of LDA+U as 
implemented in the linearized augmented planewave 
method, ■^•^ and as usual the ratio of Slater integrals 
is fixed at ^4/^2=0.688, F6/i^2=0.495. Since we 
are particularly interested in the stability of the 
atomic moment, the exchange integral J that enters 
the LDA+U method deserves attention. In atomic 
physics, and in the LDA+U method, the exchange 
integral plays two roles, ft describes the spin depen- 
dence of the Coulomb interaction, that is, the usual 
Hund's rule coupling. In addition, it carries the or- 
bital off-diagonality; with J=0 all 4/ orbitals repel 
equally by U, whereas in general the anisotropy of 
the orbitals leads to a variation^'* that is described 
by J. 

For a half filled shell for which the orbital occu- 
pations rimi = 1 and n,„j^=0 for all suborbitals m, 
the exchange effect primarily counteracts the effect 
of U , since the anisotropy of the repulsion averages 
out. As a result, using Ueff = U — J with J^ff—O is 
almost equivalent, for a perfectly half-filled shell, to 
using U and J separately as normally is done. Since 
it could be argued that Hund's first rule is treated 
adequately by the LDA exchange-correlation func- 
tional, for our calculations we have set J=0. This 
becomes approximate for the off-diagonality effects 
when the minority 4/ states begin to become occu- 
pied at high pressure. However, we have checked the 
effect at a/ao=0.8 {v = 0.5, P=60 GPa). Compar- 
ing U = 6.9 eV, J=l eV with U = 5.9 eV, J=0, 
we find the energy is exactly the same (to sub-mRy 
level) and the moment is unchanged. This result 
is in line with the Ueff, Jeff argument mentioned 
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above. Changing J from 1 eV to with U = 5.9 cV 
also leaves the energy unchanged, illustrating the 
clear unimportance of J. The J=l eV calculation 
does result in a 0.03 /is larger moment. At much 
smaller volumes, where the minority bands overlap 
the Fermi level, the changes become noticeable and 
would affect the equation of state, but only in a very 
minor way. In general, neglect of J will tend to un- 
derestimate the stability of the magnetic moment, 
which we show below already to be extremely sta- 
ble. 



B. Structure 

The observed structures of Gd follow the sequence 
hcp^Sm-type— >dhcp^dfcc -^bcm (dfcc=distortcd 
fee, which is trigonal; bcm=body-centered mon- 
oclinic). All except the bcm phase are close- 
packed arrangements, differing only in the stacking 
of hexagonal layers. The bcm phase is a lower sym- 
metry phase that suggests /-electron bonding has 
begun to contribute. 

For our purpose of studying trends relating only 
to the atomic volume, it is best to stay within a sin- 
gle crystal structure. We expect the results to refiect 
mostly local physics, depending strongly on the vol- 
ume but only weakly on the long-range periodicity. 
Therefore we have kept the simple fee structure for 
the results we present. 



III. LDA+U RESULTS 

The overall result of our study is that evolution 
of the volume and the Gd moment are predicted by 
LDA+U to be continuous under reduction of vol- 
ume, with no evidence of a volume-collapse transi- 
tion (or any other electronic phase change) in the 
region where one is observed (59 GPa), or even to 
much higher pressure. This result provides some 
support for the suggestion that the volume collapse 
is Kondo-driven, or involves in an essential way fluc- 
tuations, neither of which are accounted for in our 
approach. 

While we will usually quote volumes or the rela- 
tive volume V, it is useful to be able to convert this 
at least roughly to pressure. We provide in Fig. ^ 
the calculated equation of state, plotted as log P 
vs. V/Vo- It can be seen that the pressure is very 
roughly exponential in -V/Vo from t;=0.8 down to 
w=0.15 (2 GPa to 4 TPa). The change in slope 




FIG. 1: Log plot of the calculated pressure versus 
volume. The relatively small difference between the 
LDA+U and LDA results is evident. The relation is 
roughly exponential below V/Vo < 0.8. Current static 
diamond anvil cells will only take Gd to the V/Vo ~ 0.35 
region. 



FIG. 2; (Color online) Behavior of the calculated mo- 
ment/cell (4/ spin moment plus conduction electron 
polarization) of Gd versus reduction in volume, from 
both LDA and LDA+U methods. For the more realistic 
LDA+U method, there is very little decrease in moment 
down to V/Vo=0A5 (~110 GPa), with a rapid decline 
beginning only around V/Vo ~ 0.2 (1.5 TPa). 



around v—OA (in the vicinity of 100-200 GPa) is 
discussed below. 



A. Magnetic Moment vs. Volume 

The behavior of the total spin moment (4/ plus 
conduction) in LDA+U is compared in Fig. 13 with 
that of LDA. The general trend is similar, but the 
decrease in moment is extended to smaller volume 
by the correlations in LDA+U. Specifically, the mo- 
ment is reduced not by decrease of majority spin 
population (which would be — > /®) but rather by 
increase in the minority spin population (/'' — > 
see discussion below). Thus LDA+U enhances the 
stability of the moment by raising the unoccupied 
minority 4/ states in energy, thus reducing and de- 
laying compensation of the filled majority states. It 
has been noted elsewhere^^ that raising the minor- 
ity states is the main beneficial effect of the LDA+U 
method for Gd at ambient pressure. 

The decrease in moment is minor down to v = 0.45 
(~90-100 GPa) beyond which the decrease from 7^b 
to 6^B occurs by w = 0.2 (P - 1 TGa). Only be- 
yond this incredibly high pressure does the moment 
decrease more rapidly, as the 4/ states become band- 
like. Even in LDA this collapse does not occur until 



below [P - 300-400 GPa). With the neglect 

of fluctuations, the simplistic interpretation of the 
LDA+U results is that the Gd "bare" spin moment 
is relatively stable to ~1 TPa. 

It might be thought that, for the region of spin 
moment of 6 /is and below, where the minority oc- 
cupation is one or more, there might be an orbital 
moment of the minority system. However, at these 
volumes (see below) the minority 4/ bandwidth is 5 
eV or more, which we think makes an orbital mo- 
ment unlikely. Therefore we have not pursued this 
possibility. 



B. 4f Bandwidth 

The behavior of the 4/ states, which become 
bands, is better illustrated in Fig. O where the 
evolution of the 4/ "bands" (the 4/ projected den- 
sity of states [PDOS]) is provided graphically. At 
a/ao=0.80 (w=0.51, P « 60 GPa, where the volume 
collapse is observed) the majority PDOS is some- 
what less than 2 eV wide and still atomic-like, since 
it does not quite overlap the bottom of the conduc- 
tion band. Above this pressure range the 4/ states 
begin to overlap the conduction bands, primarily due 




FIG. 3: (Color online) View of the 4/ projected density 
of states under compression, with majority spin plotted 
upward and minority plotted downward. The curves are 
displaced for clarity, by an amount proportional to the 
reduction in lattice constant. The legend provides the 
ratio a/ao, which is decreasing from above, and from 
below, toward the middle of the figure. 



but it requires supra-TPa pressures to broaden them 
into the range of the majority 4/ states. The upturn 
in the logP vs. V curve in Fig. ^in the vicinity of 
100-200 GPa is probably due to 5p semicore overlap 
on neighboring atoms (repulsion of closed shells as 
they come into contact). The 5d bands broaden in 
the standard way under pressure, and begin to rise 
noticeably with respect to the 4/ states beyond 60 
GPa. 

The minority 4/ bands fall somewhat with respect 
to Ef as they broaden, both effects contributing to 
an increase in the minority 4/ occupation at the ex- 
pense of 5d and 6sp character. Since the majority 
4/ states remain full, the effect is that the total / 
count increases and the spin moment decreases (as 
discussed above). 

The volume dependence of the 4/ bandwidth in 
nonmagnetic Gd has been looked at previously by 
McMahan et al} They identified the intrinsic width 
Wf f from the bonding and antibonding values of the 
4/ logarithmic derivative; Wff lies midway (roughly 
halfway) between our majority and minority band- 
widths, see Fig. ^ McMahan et al. also ob- 
tained a hybridization contribution to the 4/ width; 
both of these would be included in our identified 
widths. Our widths, obtained for ferromagnetically 
ordered Gd, are difhcult to compare quantitatively 
with those of McMahan et al., because the positions 
of our minority and majority states differ by 12 eV at 
P=0, decreasing under pressure. Note that our mi- 
nority and majority widths, obtained visually from 
Fig. O differ by a factor of ~ 6 at v—l.Q, still by a 
factor of 2.5 at w=0.3, and only become equal in the 
V <Q.2 range. 



Comments on Mott Transition 



to the broadening of the conduction bandwidth. By 
a/ao=0.70 (w=0.34, P w 200 GPa) the width is at 
least 3 eV and the shape shows the effect of hy- 
bridization and formation of bands. For yet smaller 
volumes the bandwidth becomes less well defined as 
the bands mix more strongly with the conduction 
states and broaden. The minority PDOS lies in the 
midst of Gd bd bands and is considerably broader 
down to a/ao=0.70, beyond which the difference be- 
comes less noticeable. 

The position of the 4/ states relative to the semi- 
core 5p, and conduction bd states, and their evo- 
lution with volume, are pictured in Fig. ^ The 
semicore bp bands broaden to ^^10 eV by 200 GPa, 



In the simplest picture (single-band Hubbard 
model) the Mott transition is controlled by the com- 
petition of kinetic {W) and potential {U) energies, 
with the transition occurring around W U . This 
transition is normally pictured as a simultaneous 
insulator-to-metal, moment collapse, and presum- 
ably also volume collapse transition. In Gd, how- 
ever, no change in moment is observed^ across the 
volume collapse transition at 59 GPa. 

In Fig. [5l the 4/ bandwidths (majority and mi- 
nority) and the Coulomb U of McMahan et al} are 
plotted versus volume. The region [/ occurs 

around v ~0. 20-0. 25. This volume corresponds to 
a calculated pressure in the general neighborhood 
of 1 TPa, indicative of an extremely stable moment 
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FIG. 4: (Color online) Plot of band positions (lines) and 
widths (bars) of the majority and minority 4/ states, the 
semicore 5p bands, and the valence 5d bands, for ferro- 
magnetic Gd. The bar at V/Vo=0.5 (~59 GPa) marks 
the observed volume collapse transition, while the ar- 
row at 113 GPa denotes the highest pressure achieved 
so far in experiment. These results were obtained from 
LDA+U method, with U varying with volume as given 
by McMahan et aim. 



well beyond present capabilities of static pressure 
cells. This criterion however presumes a simple sin- 
gle band system, which Gd is not. 

Gunnarsson, Koch, and Martin have considered 
the Mott transition in the multiband Hubbard 
j-i-^Q(jgj^26j27j28 g^j-^j found that the additional chan- 
nels for hopping favored kinetic processes that re- 
duced the effect of the Coulomb repulsion. They 
argued that the criterion involved the inverse square 
root of the degeneracy, which can be characterized 
by an effective repulsion U* = (for / states 

the degeneracy is 2€ = 7). The Mott transition 
could then be anticipated in the range W ^ U* , for 
which U*{V) has also been included in Fig. |S1 Tak- 
ing W as the average of the majority and minority 
widths gives the crossover around Vc ^0.35 (Pc ~ 
200 GPa); taking W more realistically as the major- 
ity bandwidth gives Vc ^0.25 {Pc ^ 750 GPa). 

Another viewpoint on the "Mott transition" in the 
4/ system is that it can be identified with the 'metal- 
lization' of the 4/ bands, which might be expected to 
be where the occupied and unoccupied bands over- 
lap. These are respectively the majority and minor- 
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FIG. 5: (Color online) Plot of the 4/ bandwidths 
(both majority and minority) , together with the volume- 
dependent Coulomb repulsion U from McMahani The 
simple crossover criterion Wf « U occurs around V/Vo = 
0.20 — 0.25, corresponding roughly to a pressure of 700- 
1000 GPa. Also pictured is U* = U/V7, see text for 
discussion. 



ity bands. Significant overlap occurs only above 2 
TPa {v < 0.20) in Fig. The fact is that metalliza- 
tion (however defined for a 4/ system in the midst 
of uncorrelated itinerant conduction bands) and mo- 
ment collapse need not coincide, and the concept of 
Mott transition may need to be generalized. 



IV. SUMMARY 

In this paper we have applied the correlated band 
theory LDA-I-U method to probe the electronic and 
magnetic character of elemental Gd under pressure. 
The calculated moment decreases slowly down to 
V/Vo = 0.20 (P > 1 TPa), and only at smaller 
volumes does the moment decrease more rapidly. 
Still, no identifiable moment collapse has been ob- 
tained. Metallization, defined as overlap of unoccu- 
pied with occupied bands, also does not occur until 
the same range of volume/pressure. However, in- 
formation from studies of the multiband Hubbard 
model, and comparison of the bandwidth to U/^/N 
ratio {N~7 is the 4/ degeneracy) suggests a "Mott 
transition" might be expected in the broad vicinity 
of 500 GPa. 



The same LDA+U method, and three different 
correlated band methods have been applied to anti- 
ferromagnetic MnO. The manganese configuration is 
half-filled and fully polarized, as is Gd, with the dif- 
ference being that it is 3d and an antiferromagnetic 
insulator rather than 4/ in a background of itinerant 
bands. All methods obtained a volume collapse from 
a high-spin to low-spin configuration. Surprisingly, 
the collapse was not to nonmagnetic but rather to a 
spin-half result. 

The critical pressures for transitions suggested by 
the present study (minimum of 200 GPa, more likely 
around 750 GPa) lie well above the volume collapse 
transition that is observed at 59 GPa. At this point 
in our understanding of the 4/ shell in Gd, there 
seems to be no viable alternative to the suggestion 
by Maddox et al. that Gd provides an example of 
the Kondo volume collapse mechanism^ 



7 

V. ACKNOWLEDGMENTS 



This work has benefited greatly from a number 
of exchanges of information and ideas with A. K. 
McMahan. We have profited from many discussions 
on Gd with C. S. Yoo, B. Maddox, R. T. Scalettar, 
and A. Lazicki, and on the moment collapse question 
with M. D. Johannes and J. Kunes. We thank M. D. 
Johannes and R. T. Scalettar for a careful reading 
of the manuscript. Support from the Alexander von 
Humboldt Foundation, and the hospitality of IFW 
Dresden, during the preparation of this manuscript 
is gratefully acknowledged. This work was sup- 
ported by Department of Energy grant DE-FG03- 
01ER45876, by Strategic Science Academic Alliance 
Program grant DE-FG03-03NA00071, and by the 
DOE Computational Materials Science Network. 



^ A. K. McMahan, C. Huscroft, R. T. Scalettar, and 

E. L. Pollock, J. Computer- Aided Materials Design 

5, 131 (1998); cond-mat/9805064 
^ D. Errandonea, R. Boehler, and M. Ross, Phys. Rev. 

Lett. 85, 3444 (2000). 
^ B. Johansson, Philos. Mag. 30, 30 (1974). 

J. W. Allen and R. M. Martin, Phys. Rev. Lett. 49, 

1106 (1982). 

^ M. Lavagna, C. Lacroix, and M. Cyrot, Phys. Lett. 

90A, 210 (1982); J. Phys. F 13, 1008 (1983). 
® H. Hua, Y. K. Vohra, J. Akella, S. T. Weir, R. Ahuja, 

and B. Johansson, Rev. High Pressure Sci. Technol. 

7, 233 (1988). 

J. Akella, G. S. Smith, and A. P. Jephcoat, J. Phys. 

Chem. Solids 49, 573 (1988). 
* B. R. Maddox, A. Lazicki, C. S. Yoo, V. Iota, M. 

Chen, A. K. McMahan, M. Y. Hu, P. Chow, R. T. 

Scalettar, and W. E. Pickett, Phys. Rev. Lett. 96, 

256403 (2006). 
® J. W. Allen and R. M. Martin, Phys. Rev. Lett. 49, 

1106 (1982). 

" M. Lavanga, C. Lacroix, and M. Cyrot, Phys. Lett. 

90A, 210 (1982). 
^ K. Held, A. K. McMahan, and R. T. Scalettar, Phys. 

Rev. Lett. 87, 276404 (2001). 
^ A. K. McMahan, Phys. Rev. B 72, 115125 (2005). 

A. Svane, W. M. Temmerman, Z. Szotek, J. Laegs- 

gaard, and H. Winter, Intl. J. Quant. Chem. 77, 799 

(2000). 

^ D. Kasinathan, J. Kunes, K. Koepernik, C. V. Di- 
aconu, R. L. Martin, I. Prodan, C. E. Scuseria, N. 
Spaldin, L. Petit, T. C. Schuhhess, and W. E. Pick- 



ett, cond-mat/066543d 

B. Johansson and A. Rosengren, Phys. Rev. B 11, 
2836 (1975). 

K.Koepernik and H.Eschrig, Phys. Rev. B 59, 1743 
(1999). 

" J. P. Perdew and Y. Wang, Phys. Rev. B 45, 13244 
(1992). 

K. Schwarz, P. Blaha and G. K. H. Madsen, Comput. 
Phys. Commun. 147. 71 (2002). 

L. B. Robinson, F. Milstein, and A. Jayaraman, Phys. 
Rev. 134, A187 (1964). 

D. D. Jackson, C. Aracne-Ruddle, V. Malba, S. T. 
Weir, S. A. Catledge, and Y. K. Vohra, Rev. Sci. Instr. 
74, 2467 (2003). 
2^ D. D. Jackson, V. Malba, S. T. Wen, P. A. Baker, 
and Y. K. Vohra, Phys. Rev. B 71, 184416 (2005). 
P. Kurz, B. Bihlmayer, and S. Bliigel, J. Phys.: Con- 
dens. Matt. 14, 6353 (2002). 

A. B. Shick, A. I. Liechtenstein, and W. E. Pickett, 

Phys. Rev. B 60, 10763 (1999). 
'^'^ M. D. Johannes and W. E. Pickett, Phys. Rev. B 72, 

195116 (2005). 
25 A. B. Shick, W. E. Pickett, and C. S. Fadley, J. Appl. 

Phys. 87, 5878 (2000). 

O. Gunnarsson, E. Koch, and R. M. Martin, Phys. 

Rev. B 54, R11026 (1996). 
2^ O. Gunnarsson, E. Koch, and R. M. Martin, Phys. 

Rev. B 56, 1146 (1997). 
2* E. Koch, O. Gunnarsson, and R. M. Martin, Comp. 

Phys. Commun. 127, 137 (2000). 



